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Abstract 



We continue the analysis of the spectral curve of the normal random matrix en- 
semble, introduced in an earlier paper. Evolution of the full quantum curve is given 
in terms of compatibility equations of independent flows. The semiclassical limit of 
these flows is expressed through canonical differential forms of the spectral curve. 
We also prove that the semiclassical limit of the evolution equations is equivalent 
to Whitham hierarchy. 
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1 Introduction 



Continuing the study of the normal random matrix ensemble, presented in the 
series of papers [1,2,3,4], we explore the evolution of the spectral curve with 
respect to all the independent flows. 

In the first part of the paper, we review the basic concepts of the theory, 
relevant to many different developments, from integrability of non-linear dif- 
ferential equations to supersymmetric Yang-Mills theory [5,6,7,8,9,10,11]. 

In the second part, we derive the evolution equations and prove that their 
semiclassical limit is the universal Whitham hierarchy associated with the 
complex curve. We also notice that the Whitham hierarchy is identical to 
the set of equations which describe Laplacian growth processes - unstable 
dynamics of an interface between two immiscible phases. 



1.1 Normal matrix ensemble 



The following is a short review of normal random matrix theory. Normal ma- 
trices M commute with their hermitian conjugate: [M, M 1 "] = 0, so that both 
M and can be diagonalized simultaneously, and have complex eigenval- 
ues. The statistical weight of the normal matrix ensemble is given through a 
general potential W(M, M+) [12]: 

e^ M ' M %(M). (1) 

Here h is a parameter, and the measure of integration over normal matrices 
is induced by the flat metric on the space of all complex matrices. Using the 
standard procedure [13], angular degrees of freedom are integrated out, leading 
to the joint probability distribution of eigenvalues z±, . . . , zn, where N is size 
of the matrix: 

1 N 
j^\^n(z)\ 2 Y[^ WM) d 2 z r (2) 

Here Ajy(z) = det(z] -1 )i<jj<jv = Hi>j( z i — Zj) is the Vandermonde determi- 
nant, and 

1 . N 

= ^j/|A^)| 2 n^* ) d 2 ^ (3) 
is a normalization factor, the partition function of the matrix model. 

We consider the case where the potential W has the form 

W = -\z\ 2 + V(z) +V(z), (4) 
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where V(z) is a holomorphic function in a domain which includes the support 
of eigenvalues. We also assume that the field 



A(z) = d z V(z) 



(5) 



(a "vector potential") is a globally defined meromorphic function. 

In a proper large N limit (h — > 0, NTi fixed), the eigenvalues of matrices occupy 
a connected domain D in the complex plane, or, in general, several discon- 
nected domains. We refer to the connected components D a of the domain D 
as droplets. 

In the case of algebraic domains [4], the eigenvalues are distributed with the 
density p = — j^AW, where A = Ad z d z is the 2-D Laplace operator [14]. For 
the potential (4), the density is uniform. 

Boundary components of the droplets form a real section of a complex curve, 
defined in the following way: let us represent the boundary of the domain as a 
real curve F(x, y) = 0. If the vector potential A(z) is a meromorphic function 
(we always assume that this is the case), the function F can be chosen to be 
an irreducible polynomial. Then we rewrite it in holomorphic coordinates as 



and treat z and z as independent complex coordinates z, z. The equation 
f(z,z) = defines a complex curve. This curve is a finite-sheet covering of 
the z-plane. The single-valued function z(z) on the curve is a multivalued 
function on the z-plane. Making cuts, one can fix single-valued branches of 
this function. The boundary of the domain is a section of the curve by the 
plane where z is complex conjugate of z. 

Similarly to the hermitian matrix ensemble [5,10,11], the complex curve of 
the normal matrix ensemble is characterized by the potential W (or by the 
"vector potential" A(z)), and by a set of g + 1 integers u a (not necessarily 
positive), where g is genus of the curve. The integers are subject to the con- 
straint J2 9 a =o v a = N. If they are all positive, then they are proportional to the 
areas of the droplets of uniformly distributed eigenvalues. In this case, every 
droplet contains v a eigenvalues. 

As one varies the potential and the filling factors v a , the curve and the inter- 
face bounding the droplets evolve. Parameters of the potential (for example, 
poles and residues of the meromorphic function (5)) and filling factors are de- 
formation parameters and parameters of growth. They are coordinates in the 
moduli space of the complex curves. An infinitesimal variation of the potential 
generates correlation functions of the ensemble [14]. 




(6) 
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2 Spectral curve and wave functions 

In this section we specify the potential to be of the form (4), and introduce 
the set of biorthogonal wave functions 

i(z) = e-^)p n (4 and *„(*)= e* v «P B (z), (7) 

where the holomorphic functions Xn( z ) are ortho normal in the complex plane 
with the weight e"' 2 ' l n . Like traditional orthogonal polynomials, the biortho- 
gonal polynomials P n (and the corresponding wave functions) obey a set of 
differential equations with respect to the argument z, and recurrence relations 
with respect to the degree n. Similar equations for two-matrix models are 
discussed in numerous papers (see, e.g., [15]). 

In the basis (7), multiplication by z is represented by the L-operator (the Lax 
operator): 

L nm Xm(z) = ZXn(z) (8) 

(summation over repeated indices is implied). Obviously, L is a lower triangu- 
lar matrix with one adjacent upper diagonal, L nm = 0asm>n + l. Similarly, 
the differentiation d z is represented by an upper triangular matrix with one 
adjacent lower diagonal. Integrating by parts the matrix elements of the d z , 
one finds: 

(L^) nm X m = hd z Xn, (9) 

where U is the hermitian conjugate operator. 
The matrix elements of iJ are 

(Lt) nm = L mn = A{L nm ) + J e* w P m (z)d z P n (z)d 2 z, 

where the last term is a lower triangular matrix. The latter can be written 
through negative powers of the Lax operator. Writing d z \og P n (z) = f + 
J2k>i v k(n)z~ k , one represents in the form 

L f = A(L) + (hn)L- 1 + Y. v{k)L ~ k i ( 10 ) 

k>l 

where and (hn) are diagonal matrices with elements and (hn). The 
coefficients are determined by the condition that lower triangular matrix 
elements of A(L nm ) are cancelled. 

In order to emphasize the structure of the operator L, we write it in the basis 
of the shift operator w such that wf n = f n +\w for any sequence /„. Acting on 
the wave function, we have: 

WXn = Xn+1- 
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In the n-representation, the operators L, V< acquire the form 

L = r n w + ^Mff ) w" fc , L f = w _1 r„ + ^ w k u£\ (11) 

k>0 k>0 

Acting on Xn, we have the commutation relation ("the string equation") 

[L, L+] = ft. (12) 
This is the compatibility condition of Eqs. (8) and (9). 

Equations (11) and (12) completely determine the coefficients v!£\ r n and uffl. 
The first one connects the coefficients to the parameters of the potential. The 
second equation is used to determine how the coefficients v£\ r n and 
evolve with n. In particular, its diagonal part reads 

nft = rS-EEl«ffi/. (13) 
fc>ip=i 



2. 1 Finite dimensional reductions 



If the vector potential A(z) is a rational function, the coefficients are not 
all independent. The number of independent coefficients equals the number 
of independent parameters of the potential. For example, if the holomorphic 
part of the potential, V(z), is a polynomial of degree d, the series (11) are 
truncated at k = d — 1 . 

In this case, the semi-infinite system of linear equations (9) and the recurrence 
relations (8) can be cast in the form of a set of finite dimensional differen- 
tial equations, whose coefficients are rational functions of z, one system for 
every n > 0. The system of differential equations generalizes the Cristoffel- 
Daurboux second order differential equation valid for orthogonal polynomials. 
This fact has been observed in recent papers [16,17] for biorthogonal polynomi- 
als emerging in the hermitian two-matrix model with a polynomial potential. 
It is applicable to our case (holomorphic biorthogonal polynomials) as well. 

The semi-infinite set {xo,Xii ■ ■ ■} is then decomposed as a "bundle" of d- 
dimensional vectors 

X(n) = (Xn, Xn+l, • • • , Xn+d-lf 

(the index t means transposition, so x is a column vector). The dimension of 
the vector is the number of poles of A(z) plus one. Each vector obeys a closed 
rf-dimensional linear differential equation 

hd z x(n) = C n (z)x(n), (14) 
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where the d x d matrix C n is a "projection" of the operator onto the n-th 
(i-dimensional space. Matrix elements of the C n are rational functions of z 
having the same poles as A(z) and also a pole at the point A(oo). If A(z) is 
a polynomial, all these poles accumulate to a multiple pole at infinity. 

2.2 Spectral curve 

The semiclassical asymptotics of solutions to Eq. (14), as h — » 0, are found 
by solving the eigenvalue problem for the matrix C n (z) [18]. More precisely, 
the basic object of the the asymptotic analysis is the spectral curve [19] of the 
matrix C n , which is defined, for every integer n > 0, by the secular equation 
det(£ n (z) — z) = (here z means z ■ 1, where 1 is the unit d x d matrix). It 
is clear that the left hand side of the secular equation is a polynomial in z of 
degree d. We define the spectral curve by an equivalent equation 

f n (z,z) = a(z)det(C n (z) - z) = 0, (15) 

where the factor a(z) is added to make f n (z, z) a polynomial in z as well. The 
factor a(z) then has zeros at the points where poles of the matrix function 
C(z) are located. It does not depend on n. 

The spectral curve (15) is characterized by the antiholomorphic involution. In 
the coordinates z,z, the involution reads (z,z) \— > (z,z). This simply means 
that the secular equation det(C n (z) — z) = for the matrix C n (z) = C n (z) 
defines the same curve. Therefore, the polynomial /„ takes real values for 

z = z: 

f n (z,Z) = fn(z,z). (16) 

Points of the real section of the curve [z = z) are fixed points of the involution. 

2.3 Schwarz function 

The polynomial f n (z,z) can be factorized in two ways: 

f n (z, z) = a(z)(z - sW(z)) ...(z- Si d \z)), (17) 
where S$(z) are eigenvalues of the matrix C n (z), or 

/„(*, z) = ajz)(z - SW(z)) ...{z- Si d \z)), (18) 

where S^(z) are eigenvalues of the matrix C n (z). One may understand them 
as different branches of a multivalued function S(z) (respectively, S(z)) on 
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the plane (here we do not indicate the dependence on n, for simplicity of the 
notation). It then follows that S(z) and S(z) are mutually inverse functions: 

S(S(z)) = z. (19) 



An algebraic function with this property is called the Schwarz function. The 
equation f(z,S(z)) = defines a complex curve with an antiholomorphic 
involution. An upper bound for genus of this curve is g = (d — l) 2 , where d 
is the number of branches of the Schwarz function. The real section of this 
curve is the set of all fixed points of the involution. It consists of a number of 
contours on the plane (and possibly a number of isolated points, if the curve is 
not smooth). The structure of this set is known to be complicated. Depending 
on the coefficients of the polynomial, the number of disconnected contours in 
the real section may vary from to g + 1. If the contours divide the complex 
curve into two disconnected "halves", or sides (related by the involution), then 
the curve can be realized as the Schottky double [20,21,22] of one of these sides. 
Each side is a Riemann surface with a boundary. 

The Schwarz function on the physical sheet is a particular root, say S^(z), of 
the polynomial f n (z, z) (see (17)). It follows from (10) that this root is selected 
by the requirement that it has the same poles and residues as the potential A. 

The formal 3 semiclassical asymptote of equation (14), in the leading order in 
h, is 

X n ~etf'*#\ (20) 

Here 

dtt^ = S^dz. 

The differential dVt^ is a physical branch of the generating differential on the 
curve (see below). 

The semiclassical asymptotics was discussed in more details in [4]. To next 
order in H, it reads 

4>(z) Vdz~ ~ yjw(°°&){z) e~* (^"-Ao S(z)dz ) . (21) 
2.4 The generating differential 
The meromorphic differential 

dtt = S(z)dz (22) 
3 This formal expression ignores the Stokes phenomenon. 
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plays an important role. It is called generating differential [3,23]. On the phys- 
ical sheet it has the same poles and residues as the differential Adz: 



dtt = Adz + S-(z)dz. 



(23) 



In the following, we use the following properties of the generating differential. 

(i) The periods over a-cycles (boundaries of the droplets) are purely imaginary, 
and are integer multiples of 2ni. They compute areas of the droplets: 



The filling factors of physical droplets (belonging to the physical sheet) are 
positive. 

(ii) The real part of the integral of the differential (z — S(z))dz from some fixed 
point £ to a point on the boundary of a droplet (a point on a a-cycle) has 
the same value for all points of the boundary: 



This quantity does not depend on z, but does depend on £o unless £o is on 
the boundary. However, the difference Q — (ftp depends on the a-cycles only. 
It is equal to a b-period of the differential dVL: 



where b Q/ g is a cycle connecting a Q and a^ cycles. 
For proofs and more details, see Ref. [3]. 

Periods over b-cycles <fi a — 0o play a r °l e of chemical potentials for the fill- 
ing factors. Here denotes a chosen reference droplet. One can use chemical 
potentials to characterize evolution of the curve instead of filling factors. 



3 Evolution of the spectral curve 

3. 1 Evolution of the quantum curve 

Evolution of wave functions of the matrix ensembles with respect to a change of 
the potential is the subject of vast literature. Most of it deals with orthogonal 
polynomials of the hermitian ensemble (see, e.g., [24]). Our case is similar to 






(24) 
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the case of biorthogonal polynomials appearing in two- matrix ensembles [15]. 
We represent the evolution equations through variational derivatives S/SA(z) 
with respect to the holomorphic part of the potential. This does not imply 
any particular parametrization of the potential. The standard result reads (cf. 
[23]): 

where 

H(C) = - [log (C - L)] + - \ [log (C - L)] - \ log C, (26) 

and the notation [. . .] +j0 means the upper triangular or diagonal part of the 
matrix. This formula reflects the triangular structure of the L-operator en- 
coded by Eq.(lO) 

L f = A(L) + (Tin)L- 1 + £ v^L- k . 

k>l 

The compatibility equations 

h ik) H{z ' ) - h ih) H{z)=mz) ' H(z ' )] (27) 

determine evolution of L (i.e., the coefficients in (11)). 

Eqs. (25,27) describe an evolution with respect to deformation parameters, 
while Eq. (12) describe the growth (increasing N while keeping all harmonic 
moments fixed). 

Equations (8,25,26,27) constitute the Lax-Sato form of the Toda lattice inte- 
grable hierarchy [25] . Biorthogonal polynomials form a very particular solution 
of the hierarchy. If A(z) is a rational function, the matrices H(z) acquire a 
special structure, which is determined by poles and residues of A(z). This case 
corresponds to finite-dimensional reductions of the Toda lattice hierarchy. 



3.2 Evolution of the semiclassical curve. Whitham hierarchy 



The variation of the Schwarz function, or the generating differential (22), 
under a change of the deformation parameters is a subject of deformation 
equations. The deformation equations are already built in the decomposition 
of the Schwarz function (23). 

In order to write them, we need holomorphic and meromorphic differentials 
canonically normalized with respect to the a-cycles; they are defined in the 
following: 
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(i) Holomorphic differentials (Abelian differentials of the first kind) dW a (z), 

(ii) Meromorphic differentials (Abelian differentials of the third kind) dW^'® (z) 
having poles at z = ( and its mirror on the back side with the residues 
±1, normalized so that § a dW^'^(z)dz = 0. The meromorphic differential 
dW(°°'°°\z) plays a special role. 

The variation of dft = Sdz with respect to t — UN gives the meromorphic 
differential dW {co '™\z): 

d t S{z)dz = -dW {o °'™\z). (28) 

This follows from the fact that the singular part of the Schwarz function (A(z) 
in (23)) does not depend on t, and that S-(z) — > t/z at infinity. 

Variations with respect to the filling factors = %v a at a fixed t and a fixed 
potential affects only the S(z) — t/z part of the Schwarz function. It gives 
some g- holomorphic differentials. We call them dW a (z): 

d t ( a )S(z)dz = —dW a (z). 

These differentials are canonically normalized: f dWp = 5 a /3, so they are 
unique. 

Finally, a variation of dQ = Sdz with respect to A(Q, at fixed filling factors, 
produces a meromorphic differential which has simple poles at z = (, oo and 
their mirrors on the back side. The poles come from the variation of the first 
term of (23). Their residues are ±1 and the sum of the residues is zero on each 
side of the double. Since all filling factors are fixed, this differential has zero 
a-periods and we get 

r 

——S(z)dz = -dW {a \z) + dW {oo ^\z), on the front side. (29) 

This equation simply reflects the fact that all singularities of the Schwarz 
function are the singularities of A(z). 

Compatibility equations constitute exchange relations among the differentials. 
They have the form of the Whitham hierarchy [26] : 

v(C)w^(z) = v(^(«'(4 

< V(()W(°°'™\z) = d t W^(z), (30) 

V(()W a (z) = d t(a) W^(z), 

V 

where V(z) = d t + ^ + 
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Now we can check that the Whitham equations (30) are the semiclassical limit 
of the quantum deformation equations (27). Substituting the semiclassical 
form of the wave function (21) into (25), we see that the classical limit of the 
flows are primitive functions of the meromorphic differentials dW^'^(z): 

Similarly, the classical limit of the shift operator w defined in Sec. 2 is equal 
to pW 00 ' 00 )^), and the second equation of (30) is a classical limit of the com- 
patibility conditions. 

The deformation equations become more illustrative for algebraic domains 
[3] and especially significant for simply-connected algebraic domains [1]. In 
these cases all the differentials are expressed through the Green function for 
the exterior Dirichlet boundary problem. Differentials of the third kind are 
expressed through the Green function on the physical sheet, 

dW {a \z) =2d z G((,z)dz, 

while the holomorphic differentials are related to harmonic measures of the 
droplets: 

dW a (z) = d z uj a {z)dz. 



We recall that the Green function G((, z) is a symmetric and harmonic func- 
tion everywhere in the exterior of D except at z — £, where it has a logarithmic 
singularity, G((,z) ~ log |C — z\. The Green function vanishes if any of the 
arguments belongs to the boundary. The harmonic measure of a droplet D a 
is a harmonic function in the exterior of D such that it is equal to 1 on the 
boundary of D a and vanishes at boundaries of all other droplets. 

If there is only one physical droplet, and no virtual droplets, the Green function 
is given by a conformal map w(z) of the exterior of the droplet to the exterior of 
the unit disc: G(C, z) = log ™(0-™( z ) _ j n ^is case ^ deformation equations 
read 

V(()S(z) = -2d z G((,z), d t S{z) = d z \ogw(z). (31) 

The second equation in (31) describes the evolution of the conformal map 
when the area of the domain increases, while all the harmonic moments stay 
constant. This equation, in different equivalent forms, has been known in the 
theory of Laplacian growth for a long time [27]. The first equation (31) ap- 
peared in [1,23,28]; for its generalization for multiply-connected algebraic do- 
mains, see [29] . A generalization of deformation and growth equations for the 
Laplacian growth of multiply-connected domains was reported in Ref. [3]. It 
appeares that they constitute the same Whitham hierarchy as the semiclassical 
limit of the evolution of the normal matrix ensemble (30). 
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4 Appendix. Laplacian Growth 



Laplacian growth referres to growth of a planar domain, whose boundary 
propagates with a velocity proportional to the gradient of a harmonic field. 
The Hele-Shaw problem is a typical example. It describes the dynamics of a 
2D system of two immiscible fluids (such as oil and water). The fluids are 
confined between two horizontal glass plates, separated by a small distance, 
so that the problem is essentially two-dimensional. 

The motion of the interface follows from the Navier-Stokes equation being 
specified for a cell. It reads: velocity field of oil is proportional to gradient of 
pressure. The liquids are assumed to be incompressible, so that the pressure 
is a harmonic function. 

In Refrs. [1,2,3], we identified the Laplacian growth with the Whitham hier- 
archy. It appeares to be identical to equations (30) describing the growth and 
deformation of the support of eigenvalues of the normal matrix ensemble. 
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